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Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of 
three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus 
domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus 
porosus) were investigated using circular dichroism, incoherent elastic neutron scattering 
and coarse-grained Brownian dynamics simulations. The experimental results from Hb sol- 
utions revealed a direct correlation between protein resilience, melting temperature and 
average body temperature of the different species on the 0.1 ns time scale. Molecular forces 
appeared to be adapted to permit conformational fluctuations with a root mean square dis- 
placement close to 1.2 A at the corresponding average body temperature of the endotherms. 
Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit 
over the entire temperature range in which the animal lives. In fully hydrated powder samples 
of human and chicken, Hb mean square displacements and effective force constants on the 
1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, 
in contrast to the solution case. A complementary result of the study, therefore, is that 
one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the 
pico- to nanosecond time scale which might be relevant for biological function. Coarse- 
grained Brownian dynamics simulations permitted to explore residue-specific effects. They 
indicated that temperature sensing of human and chicken Hb occurs mainly at residues 
lining internal cavities in the (3-subunits. 
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circular dichroism 



1. INTRODUCTION 

Haemoglobin (Hb) is the major macromolecular com- 
ponent of red blood cells (RBCs). The protein forms a 
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tetramer under physiological conditions of the two so- 
called a- and two (3-chains. Each subunit acts as a scaf- 
fold for one haem-group, which can bind oxygen 
reversibly. One of the biological functions of RBCs is 
to transport oxygen from the lung to the tissues. Only 
higher vertebrates produce these highly specialized 
cells for oxygen transport in their bodies. Higher 
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vertebrates can be subdivided into two sub-classes: 
endotherms — such as mammals or birds — which maintain 
a constant body temperature mainly through regulation 
of their internal metabolism, and ectotherms — such as 
reptiles or amphibians — whose body temperature is 
regulated by the environment. 

The mechanical properties of single human RBCs 
have been studied by micropipette aspiration exper- 
iments as a function of temperature [1]. A passage 
transition occurring within a very narrow temperature 
range has been found in these experiments. Below a 
temperature of 36.4 + 0.4°C, all aspirated RBCs bloc- 
ked the pipette tip (diameter of the pipette approx. 
1.3 jxm) and were just compressed in volume. The Hb 
concentration within the RBCs, in this case, reached 
very high values up to approximately 0.5 g ml -1 , 
whereas under physiological conditions, the cellular 
Hb concentration is approximately 0.33 g ml -1 [1]. 
Above the passage transition temperature, all aspirated 
cells entered the micropipette tip easily without any 
apparent resistance [1]. The passage transition tempera- 
ture was surprisingly very close to human body 
temperature of approximately 37° C, which suggested 
a physiological relevance. Consecutively, the viscosity 
of highly concentrated Hb solutions from 0.33 to 
0.5 g ml -1 was measured as a function of temperature 
[1]. A drop in the viscosity of the highly concentrated 
Hb solutions was found close to the passage transition 
temperature reminiscent of a colloidal gel-to-fluid tran- 
sition. Hence, the passage transition of the RBCs was 
attributed to a more fluid Hb solution in the cells 
above the body temperature. Colloid osmotic pressure 
measurements on whole RBCs showed Hb aggregation 
starting at body temperature [2]. Small angle neutron 
scattering experiments revealed the presence of a 
large-scale superstructure with a mass fractal dimension 
of 2.67 in a concentrated Hb solution of 290 mg ml -1 , 
which might provide the basis for body temperature- 
related Hb aggregation [3]. Further experiments have 
been performed to elucidate the molecular origin of 
the observed effects. Changes of secondary structure 
content of human Hb around the passage transition 
temperature were monitored by circular dichroism 
(CD), which revealed a significant loss of a-helical con- 
tent of human Hb at 37.2 + 0.6°C being again close 
both to body temperature and to the observed passage 
transition temperature [4]. The partial loss of secon- 
dary structure close to body temperature has been 
introduced as 'partial unfolding' in the literature. Fur- 
thermore, it has been investigated whether Hb 
molecules from endothermic vertebrates would show 
similar partial unfolding transitions close to their 
specific body temperatures [5]. Interestingly, this was 
indeed the case. The partial unfolding temperature of 
Hb from a large variety of species was directly corre- 
lated with the body temperature of the animal, 
ranging from 34°C for the duckbilled platypus, which 
was the animal with the lowest body temperature in 
that study, to 42°C for a bird, the spotted nutcracker, 
which had the highest body temperature in the study 
[5]. Incoherent quasi-elastic neutron scattering exper- 
iments on human RBCs and on concentrated human 
Hb solutions showed that the partially unfolded state 



is more flexible and softer than the folded state below 
body temperature [6,7]. In that sense, Hb appears to 
act as a molecular thermometer and to sense body 
temperature [2] , which would be of outstanding interest 
and importance for biological, bioengineering and bio- 
medical RBC research. As both the secondary 
structure content and the three-dimensional structure 
of Hb are highly similar [5] , the mechanism underlying 
temperature sensing needs to be explored in the specific 
physico-chemical properties of the different amino acids 
and in their interactions [2]. One hypothesis would be 
that stronger internal forces in the structure of Hb 
would be responsible for the observed correlation 
between partial unfolding transitions and body temp- 
eratures. Thus, stronger internal forces might stabilize 
the structure of Hb and would be tuned very finely to 
respond to small differences in body temperature. How- 
ever, to our knowledge, no studies are reported in the 
literature that investigated the relation between protein 
thermal stability, internal forces and conformational 
fluctuations in various Hb. 

The aim of our study was to investigate whether 
internal dynamics and thermodynamic parameters of 
Hb from different species are adapted to body tempera- 
ture, and to identify the thermo- sensitive amino acids 
on a residue- specific level. Neutron spectroscopy is a 
well-suited experimental technique to measure average 
fluctuations of proteins in the pico- to nanosecond time 
scale and in absolute units of a few angstroms [8-11]. 
Although significant differences were found in the internal 
dynamics of Hb in highly concentrated solution, the mean 
square displacements (MSD) of hydrated human and 
chicken Hb powder samples showed no differences over 
the whole temperature range from 10 to 300 K demon- 
strating that not all biologically relevant motions of Hb 
are activated with only one hydration layer. Coarse- 
grained Brownian dynamics simulations on the other 
hand are a complementary and very informative tool as 
they provide residue-resolved information on protein flexi- 
bility [12-14] . In order to obtain a detailed understanding 
of the experimental results and to explore residue-specific 
effects, coarse-grained Brownian dynamics simulations 
were performed on human and chicken Hb. 



2. MATERIAL AND METHODS 
2.1. Sample preparation 

Human Hb was purchased from Sigma-Aldrich. Blood 
samples of platypus, chicken and salt water crocodile 
were obtained from live animals. RBCs were harvested 
by centrifugation and washed several times with isotonic 
buffer. The platypus RBCs were then lysed by adding 
distilled water and centrifuged at 20 000 relative cen- 
trifugal force (rcf). Crocodile and chicken RBCs are 
nucleated and were lysed with a lysis buffer (0.5% Tris- 
HC1 at pH 7.5) to avoid breaking the nuclear membrane. 
Nuclei were separated by centrifugation at 2800 rcf and 
the supernatant was spun at 20 000 rcf. The Hb solutions 
were dialysed against distilled water before lyophiliza- 
tion. For the neutron scattering experiments, the Hb 
powders were dissolved in 99.9 per cent D 2 0 to exchange 
the labile H atoms and freeze dried afterwards. 
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For the experiments on the IN 13 neutron spectro- 
meter (see §2.2), highly concentrated Hb solutions were 
prepared. The dry Hb powders were placed in flat alu- 
minium sample holders of beam path 0.20 or 0.30 mm. 
D 2 0 buffer (0.1 M KC1, 61.3 mM K 2 HP0 4 , 5.33 mM 
KH2PO4, pD 7.4) was then added to a level of 1.1 g 
D 2 0/1 g Hb. The following amounts of Hb were used 
for the experiments on INI 3: 23 and 75 mg platypus 
Hb, 118 mg crocodile Hb, 250 mg chicken and human 
Hb. Hydrated Hb powder samples were prepared for the 
experiments on the SPHERES neutron spectrometer 
(see §2.2). Human and chicken Hb were placed in flat alu- 
minium sample holders and dried in a desiccator over 
silica gel until no further loss of weight was observed. 
The dried Hb powders were rehydrated to a level of 
0.4 g D 2 0/1 g Hb in a saturated D 2 0 atmosphere. The 
hydration level was determined by weighing the samples. 
Around 300 mg of human and chicken Hb were used for 
the experiments on SPHERES. The beam path of the 
sample holders was 0.30 mm. 

Hb samples for CD spectroscopy were prepared at a 
concentration of 0.3mgml _1 in 20 mM K 2 HP0 4 / 
KH 2 P0 4 buffer at pH 7.0. 



2.2. Neutron scattering experiments 

Neutron scattering was measured on the thermal neutron 
backscattering spectrometer IN13 at the Institut Laue- 
Langevin in Grenoble, France [15], and on the cold 
neutron backscattering spectrometer SPHERES at the 
FRM II in Garching, Germany [16]. The relation between 
the modulus of the scattering vector q, the incident neutron 
wavelength A and the scattering angle 26 is q=AiT/X 
sin( 0) for neutrons without energy transfer. IN 13 is charac- 
terized by an energy resolution of AE = 8 jxeV (FWHM) , 
and a scattering vector range of 0.2 < q < 4.9 A -1 . 
SPHERES is characterized by a very high-energy resol- 
ution of AE=0.66 |xeV (FWHM) and a scattering 
vector range of 0.17 < q< 1.84 A -1 . IN13 is sensitive 
to macromolecular motions in the space-time window 
of a few angstroms in approximately 0.1 ns, whereas 
SPHERES detects slower fluctuations with amplitudes of 
a few angstroms in approximately a few nanoseconds. 

The neutron detectors of IN13 were calibrated with a 
vanadium measurement at 280 K. On SPHERES, the 
Hb powder samples were cooled to 10 K, and the low- 
temperature measurements were used for calibration. 
The incoherent contribution of D 2 0 in the concentrated 
Hb solution samples was estimated to be less than 4 per 
cent of the total incoherent scattering [7] and is negligible 
for the hydrated powder samples. On IN13, the tempera- 
ture was varied stepwise, while on SPHERES, the 
temperature was ramped up continuously from 10 to 
300 K within 1 day, and measured data were binned 
into temperature steps of 5 K. Measured data of the 
empty sample holder were subtracted from the samples. 
No multiple scattering corrections were performed as 
the transmissions of all samples were above 0.9. 

Mean square displacements (u 2 ) were calculated 
within the Gaussian approximation according to 



where /( q) represents the measured elastic intensity [9] . 
The (u 2 ) were determined from the initial slope of 
the measured elastic intensity in the q 2 - range of 1.6 < 
/<3.5 A" 2 for IN13, and 0.36 < q 2 < 3.10 A" 2 for 
SPHERES, where the initial slope of the curves was 
found to be linear. The Gaussian approximation of 
protein dynamics [17] is mathematically identical to 
the Guinier approximation of small-angle scattering 
[18,19]. The Gaussian approximation in the strictest 
sense holds up to (u 2 ) • q 2 ^2 for confined motions of 
all shapes [9]. However, protein dynamics is intrinsically 
anisotropic with motions of ellipsoidal shape and axial 
rations of 1:1:2 [20] . The Guinier approximation 
holds to larger values for such ellipsoidal shapes [19]. 
For protein dynamics, the Gaussian approximation 
was found to be valid for (u 2 ) ■ q 2 -S [21,22]. The (u 2 ) 
values measured with IN 13 in our study are in the 
range up to q 2 • (u 2 ) < 6.8. For SPHERES, they are 
below q 2 ■ (u 2 ) < 3.5. 

Effective mean force constants (kf) that describe the 
mean resilience of the proteins were obtained from 
the slope of (u 2 ) versus temperature T according to 



(kf): 



0.00276 
d(u 2 )/dT' 



(2.2) 



(u 2 ) 



-6- A[ln/(g)] 



(2.1) 



The units are chosen such that (kf) is in N m 1 when 
(u 2 ) is in A 2 and T is in K [9,23]. 

2.3. Circular dichroism experiments and 
thermodynamic analysis 

A J-810 spectropolarimeter (JASCO Corp., Tokyo, 
Japan) equipped with a temperature-controlled cuvette 
holder was used to record the CD spectra . The tempera- 
ture-dependent loss of a- helical structure was determined 
from the CD signal at a constant wavelength of 222 nm. 
The samples were measured in 1 mm optical path 
quartz cuvettes under constant N 2 flow at a heating rate 
of 1°C min -1 . The raw data were corrected for pre- 
and post-transitional slopes. A van't Hoff analysis was 
used to determine the thermodynamic values of AH = 

^unfolded — ^folded an( ^ AS = S^nfoi^ed — Sfoided between 

the unfolded and folded state at the melting tempera- 
ture T m . The model assumes a two-state equilibrium 
between the unfolded and the folded state. In the unfold- 
ing transition region, the difference in free-energy AG 
was calculated using AG{T) = —RT ln(/u//r) with the 
folded fraction / F , the unfolded fraction fj= 1 — f F and 
the molar gas constant R. The enthalpy and entropy 
differences at T m were determined according to 
AG( T) = AH - T ■ AS in the linear regime. 

2.4- Sequence alignments 

Protein amino acid sequences were obtained from the 
NCBI gene and protein databases (http://www.ncbi. 
nlm.nih.gov/). Multiple sequence alignments were 
performed using the Clustal W software [24]. The 
a-helical content of the subunits was calculated using 
the Jpred3 secondary structure prediction server [25]. 
Entry names and accession numbers of the sequences 
of the alpha- and bet a- units are: human Hb, 
HBA_HUMAN P69905 and HBB_HUMAN P68871; 
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platypus Hb, HBA_ORNAN P01979 and HBB_OR- 
NAN P02111; chicken Hb, HBA_CHICK P01994 and 
HBB_CHICK P02112; crocodile Hb, HBA_CRONI 
P01998 and HBB_CRONI P02129. 

2.5. Computer simulations 

Coarse-grained Brownian dynamics simulations were 
carried out on the biological assemblies of the human 
and chicken Hb deoxy crystal structures (pdb access 
codes 2DN2 and 1HBR) using the ProPHet (probing 
protein heterogeneity) program [13,26,27]. The simu- 
lations used a coarse-grained protein model, in which 
each amino acid is represented by one pseudoatom 
located at the C a position, and either one or two (for 
larger residues) pseudoatoms replacing the side chain 
(with the exception of Gly) [28]. Interactions between 
the pseudoatoms were treated according to the stan- 
dard elastic network model [29]. The elastic network 
model is a simplification of the heterogeneity of internal 
protein forces as all pseudoatoms lying closer than 9 A 
are joined with quadratic springs having the same 
force constant of 0.6 kcal mol -1 A -2 = 0.42 N m -1 . 
Springs are assumed to be relaxed in the crystallo- 
graphic structure of the protein. Following earlier 
studies, which showed that the haem group had little 
influence on calculated force constants [13,27], we 
chose not to include the prosthetic group in the protein 
representation. The simulations used an implicit solvent 
representation via the diffusion and random displace- 
ment terms in the equation of motion [30], and 
hydrodynamic interactions were included through the 
diffusion tensor [31]. Further details regarding the 
simulation procedure can be found in Sacquin-Mora & 
Lavery [13] and Sacquin-Mora et al. [27]. 

The Brownian dynamics simulations were performed 
with 200 000 steps at an interval of 10 fs and at a temp- 
erature of 300 K. Effective force constants for displacing 
each particle i were calculated as 

k= 3fer 2 , (2.3) 

where the brackets indicate the average taken over the 
whole simulation, fee is the Boltzmann constant and di 
is the average distance of particle % from the other particles 
j in the protein, excluding the pseudoatoms, which belong 
to the same residue m to which also particle i belongs. The 
distances between the C a pseudoatom of residue m and 
the C a pseudoatoms of the adjacent residues m + 1 and 
m — 1 are not included in the average. The force constant 
associated with each residue m is taken to be the average 
of the force constants calculated according to equation 
(2.3) for each of the pseudoatoms % forming this residue. 
Within this framework, the mechanical properties of the 
protein are described at the residue level by its 'rigidity 
profile', that is, by the ordered sequence of the force 
constants calculated for each residue. 



3. RESULTS AND DISCUSSION 

Mechanisms of thermo-adaptation have been studied 
mainly for proteins from bacteria and archaea, which 



live at extremes of temperatures and are specially 
adapted to their optimal living conditions. Adaptation 
to high temperature can be achieved in various ways, 
for example, by stronger hydrogen bonding, clusters of 
charged amino acids or reduced internal cavities [32]. 
Hb on the other hand is a protein from higher ver- 
tebrates and it is not necessarily the case that 
mechanisms found in extremophiles determine thermo- 
sensation in Hb. 

3.1. Multiple sequence alignments and secondary 
structure prediction 

The structures of human and chicken Hb have been 
solved by X-ray crystallography [33,34]. Both proteins 
are highly similar in their structures. The crystal struc- 
tures of salt water crocodile and platypus Hb have not 
yet been determined. We performed multiple sequence 
alignments using the Clustal W software [24] and pre- 
dicted the a-helical regions from the aligned amino acid 
sequences using the JPred3 secondary structure predic- 
tion server [25]. The amino acid sequence of salt water 
crocodile is not published. As a compromise, we used 
the amino acid sequence of nile crocodile, which is a 
very close relative [35]. Both crocodile proteins are 
expected to be similar in their sequences and in their 
secondary structures. The aligned sequences and the 
predicted a-helical regions of all proteins are shown in 
figure 1. Seven a-helical regions are predicted both for 
the a- and (3-subunits of all Hb proteins. Note that 
the crystal structures show eight helices. The length 
and location of the predicted regions are nearly identi- 
cal for all proteins. The sequence alignments 
demonstrate that all proteins share a large amount of 
similar and identical amino acids. However, there are 
also a considerable number of non-similar amino acids 
distributed in the sequences. Regions with a larger 
number of non-similar amino acids are located before 
and in the first predicted helix of the a- and (3-subunits, 
in the small loop between the predicted helices F and G 
of the a-subunit and in the large loop between the pre- 
dicted helices C and D of the p-subunit. Although, a 
high similarity of the secondary structure of all proteins 
was found, we cannot infer from the sequence align- 
ments which residues are responsible for differences 
occurring in the thermodynamic properties and in the 
conformational fluctuations. Protein thermal motions 
and thermodynamics depend on the local heterogeneous 
chemical environment and also on the interactions of 
different amino acids which might be separated in the 
sequence, but which are actually in close contact in 
the structure. 

3.2. Thermodynamic analysis 

The thermal stability of Hb was studied by measuring 
the unfolding transitions using CD. We assumed a 
direct proportionality between the folded fraction of 
the protein and the recorded ellipticities. The obtained 
unfolding curves are shown in figure 2 a, and the derived 
free-energy difference A G = G U nfoided ~~ Glided between 
the unfolded and the folded states in the transition 
region is given in figure 2b. The enthalpy and entropy 
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(a) 

HbA_human 
HbAplatypus 
HbA_crocodile 
HbA chicken 



HbA human 
HbA_platypus 
HbAcrocodile 
HbA chicken 



MVLSPADKTNVKAAWGKVGAHAGEYGAEALERMFLS;FPTTKTYFPHFDLSHGS@QVKGHG 
-MLTDAEKKEVTALWG^AGH^EYGAEALERLF^AFPTTKTYFSHFDLSHGSraiKAHG 
-VLSSDDKCNVKAVWS^GHIEEYGAEALERMFWPQTKIYFPHFDLSHGSAQIIRAHG 
-MLTAEDKKLIQQAWEKA(ASH3EEFGAEALTRMFT^YPQj|KTYpP|FDLSPGSDQVRGHG 



.* #* ** ***** ** *. 



KKVADALTNAVA|VDDMPNALSALSDiHA|KLRVDPVNFKLLSHCLLVTLAAlLPAEFTP 
KKVADALSTAAGHFDDMDSALSALSDLHAHKLRVDPVNFKLLAHCILVVLARHCPGEFTP 
KKVFAALHEAVN|IDDLPGALCRLSELHAHSLRVDPVNFKFLA|QCVLVVVAIHHPGSLTP 
KKVLGALGNAVKNVDNLSQAMAELSNLH^YNLRVDPVN[FKL|LSQCIQVVLAVHMGKDYTP 



*** ** 



HbAhuman 
HbAplatypus 
HbA_crocodile 
HbA_chicken 

(b) 

HbBhuman 
HbB_platypus 
HbBcrocodile 
HbB chicken 



HbBhuman 
HbB_platypus 
HbBcrocodile 
HbB chicken 



AVHASLDKFLASVSTVLTSKYR 
SAHAAMDKFLSKVATVLTjSKYR 
EVHASLDKFLCAVSSVLTSKYR 
EVHAAFDKFLSAVSAVLAE|KYR 



.**. *** 



mvhltpeeksavtalwgkvnvd^ggealgrlp/v|ypw frqrf| fesfgdlstpdavmgn^ 
-vhlsggeksavtnl^gkvnineijggealgrllvvtypwtqr^feafgdlssagavmgnpk 
-asfdphekqligdlwhkvdvahlcggealsrmlivypwkrr^fenfgdisnaqaimhnek 
mvhwtaeekqlitglwgkvnvae|c^ealarllivyp ^tqrf| fasfgnlssptailgnpm 



vkahgkkvlgafsdglahldnlkgtfatlselhci^klhvdpenfrllgnvlvcvlahhjfg 
vkahgakvltsfgdalklldd^kgtfaklselhcdklhvdpenfnrlgnvlivvla^fs 
vqahgkkvlasfgeavchldgirahfanlsklhcp5lhvdpen]fkllgdiiiivlaah)yp 
vrahgkkvltsfgdavkPldnikntfsqlselhcdklhvdpe^f]rllgdiliivlaa(kfs 



* . *** *** 



**.***.********* 



HbBhuman 
HbBplatypus 
HbB_crocodile 
HbB chicken 



KEFTPPVQAAYQKVVAGVANALAlHlKYH 
KDFSPEVQAAWQKLVSGVAHALGHKYH 
KDFGLECHAAYQKLVRQVAAALAAEYH 

KDFTPECQAAWQKLVRVVAHALAgKYH 

*.* .**.**.* ** ** . ** 



Figure 1. Amino acid sequence alignment of the (a) a- and (b) P-subunits of human, platypus, crocodile and chicken Hb. 
The amino acid similarity between all proteins is given as follows: asterisks indicate identical amino acids, colons indicate similar 
amino acids and empty space represents absence or low degree of similarity. The grey background indicates predicted a-helical 
contents of the subunits. Black frames label the parts of predicted helices where the prediction confidence was less than 50%. 
(Online version in colour.) 
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Figure 2. (a) Unfolding transitions of the investigated Hb 
samples, (b) Difference in free-energy AC between the unfolded 
and the folded state of the proteins in the transition region. 
Squares, platypus Hb; circles, human Hb; triangles, crocodile 
Hb; inverted triangles, chicken Hb. (Online version in colour.) 



differences AH and AS at the melting temperatures T m 
were determined from AG. The thermodynamic par- 
ameters T m , AH and AS are summarized in table 1. 
The T m of chicken and human Hb are in good agree- 
ment with melting temperatures of avian Hb (goose 
68.2°C, chicken 67.6°C and turkey 67.8°C) and 
human Hb (64°C) as reported by Ajloo et al. [36]. To 
our knowledge, no melting temperatures of Hb of 
salt water crocodile or platypus are reported in the 
literature. Zerlin et al. [5] reported average body temp- 
eratures T h of healthy platypus, human and chicken of 
33.0 ± 1°C, 36.6 ± 0.2°C and 41.0 ± 0.5°C, respect- 
ively. For comparison, the body temperatures are also 
given in table 1. The errors of the body temperatures 
refer to statistical variations during repeated measure- 
ments of healthy animals under the same conditions. 
Obviously, the body temperature of an animal depends 
on its physical activity, on the place of the body where 
the measurement was made and even varies in the cir- 
cadian cycle. For example, the body temperature of 
freely living platypuses was found to vary within 
+ 3°C around the average value [37]. Body temperature 
can also vary between different animals of the same 
species. A statistical analysis of literature data was 
reported 34.4-37.8°C as a typical body temperature 
range of humans [38]. From a cautious point of view, 
the body temperature range of the different species 
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Table 1. Average body temperatures T b , thermodynamic parameters T m , AH and AS of the different Hb samples. The 


enthalpy and entropy differences between the unfolded and the folded states 


were determined at i m . ihe 


errors of T m are 


around 0.5%, the errors or AH and Ab are around 3%. 


The statistical errors 


ol i b ol the endotherms are 


below 3%, body 


temperature ranges of up to + 3°C are reported. 








platypus Hb 


human Hb 


chicken Hb 


crocodile Hb 


T b (°C) 33 


36.6 


41 


25-34 


T m (°C) 60.5 


63.8 


68.7 


68.3 


Ai^kJmol" 1 ) at T m 177 


320 


471 


162 


AS (kJ mol" 1 K" 1 ) at T m 0.529 


0.950 


1.379 


0.476 



should be considered more as + 3°C around the average 
value. The measured melting temperatures T m of the 
three endotherms show a direct correlation with 
the reported average body temperatures (see figure 4). 
A similar correlation between body and environmental 
living temperatures was reported in the literature for 
eye lens crystallins of different vertebrates, ranging 
from an average living temperature of — 1.9°C and a 
melting temperature of T m = 47° C of the Antarctic tel- 
eost fish to an average living temperature of 47.0°C and 
a melting temperature of T m = 67.4°C of the desert 
iguana lizard [39] . Crocodiles are in that sense a special 
case. As ectotherms, they do not metabolically main- 
tain a constant body temperature. Instead, they 
regulate their body temperature by moving between 
water and land. In this way, they can keep their body 
temperature in a relatively narrow range by their behav- 
iour. Grigg et al. [40] reported a body temperature 
range of salt water crocodiles between 25 °C and 34° C. 
The melting temperature of crocodile Hb was found to 
be rather high, within the error bars identical to chicken 
Hb, although the physiological body temperature range 
is situated at comparatively low temperatures. 

Further information on the mechanisms of thermal 
stabilization can be obtained from the thermodynamic 
parameters AH and AS at T m . Concerning Hb from 
the endotherms, higher thermal stability is correlated 
with an increase of both AH and AS at T m . A larger 
value of AiJ(T m ) shifts the free-energy difference AG 
to larger values and increases thermal stability. An 
increase of AS( T m ) at the same time results in a steeper 
slope of AG, which results in a reduction of the melting 
temperature (figure 2b). As a consequence, the enthal- 
pic component seems to dominate the thermodynamic 
stabilization of Hb from the endotherms. A different 
mechanism of thermal stabilization can be identified 
for crocodile Hb: Crocodile Hb has a high- melting temp- 
erature similar to that of chicken Hb, but the enthalpy 
and entropy differences at the melting temperature of 
crocodile Hb are nearly three times smaller than those 
of chicken Hb (compare table 1). Intriguingly, the 
values of AH(T m ) and AS(T m ) of crocodile Hb are 
even smaller than those of platypus Hb, which has the 
smallest melting temperature of the investigated pro- 
teins. The small value of Ai7(T m ) of crocodile Hb 
destabilizes the folded state of the protein. The enthal- 
pic destabilization is compensated by a small AS(T rn ) 
value, which stabilizes the folded state. The smaller 
entropy difference A S at T m of crocodile Hb compared 
with that of platypus Hb results in a flatter slope of 
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Figure 3. Mean square displacements (u 2 ) of the highly con- 
centrated Hb solutions as a function of temperature. Black 
circles, platypus Hb; red triangles, human Hb; green inverted 
triangles, crocodile Hb; blue diamonds, chicken Hb. 



AG in the transition region, which shifts the T m of cro- 
codile Hb to a higher value. Indeed, the combination of 
enthalpic and entropic terms to protein stabilization is 
subtle and carefully counterbalanced for crocodile Hb 
to reach a high- melting temperature. 



3.3. Incoherent elastic neutron scattering 

The Hb samples as highly concentrated solutions (1.1 g 
D 2 0/1 g dry Hb) were measured using the thermal neu- 
tron backscattering spectrometer INI 3. The instrument 
is sensitive to macromolecular conformational motions 
on the time scale of a few 0.1 ns and with amplitudes of 
a few angstroms. The average MSD, (u 2 ), were determined 
within the Gaussian approximation from the recor- 
ded intensities. Representative measured data of all 
samples are given in electronic supplementary material, 
figure SI. The determined MSD of all samples are 
shown in figure 3. Only a comparatively small amount 
of platypus Hb was available for the experiments, which 
explains the large errors. The measurement of platypus 
Hb was repeated with a sample prepared from a different 
animal as a control. The MSD of both samples overlap 
within the error bars, and the (u 2 ) were combined for 
the analysis. Effective force constants (kf) and the MSD 
at the reported body temperatures [5] were determined 
from linear fits to the (u 2 ). The obtained values of (kf) 
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Table 2. Force constants and root mean square displacements obtained from neutron scattering. Effective force constants (kf) 
and RMSDs at T b of the concentrated Hb solution samples (IN13 data), and effective force constants (ki), (A^), (k£) of 
the hydrated powder samples (SPHERES data). The values of (hi), (ki) and (k%) were calculated in the temperature regions 
10-100 K, 100-200 K and 252-292 K, respectively. 



platypus Hb 



human Hb 



chicken Hb 



crocodile Hb 



IN13 

(kf) (N m" 1 ) 
RMSD at T b 
SPHERES 
(h) (N m" 1 ) 
(hi) (N m" 1 ) 
(k{) (N m" 1 ) 



(A) 



0.11 ± 0.01 
1.39 ± 0.11 



0.15 + 0.02 
1.22 + 0.03 



6.8; 
1.0; 
0.21 ; 



; 1.2 
0.1 

0.02 



0.23 ± 0.03 
1.23 ± 0.04 

5.6 + 0.8 
1.1 + 0.1 
0.24 + 0.02 



0.18 ± 0.02 




64 66 

Figure 4. Melting temperatures T m , effective force constants (kf) 
and body temperatures T b of the endotherms. The dashed line 
indicates a direct correlation between the results. Filled squares, 
effective force constants; open circles, body temperature. 



and the val ues o f the root mean square displacements 
(RMSD), yj (u 2 ), at the body temperatures of the three 
endotherms are reported in table 2. 

Tehei et al. [21] studied average macro molecular 
dynamics in psychrophilic, mesophilic, thermophilic 
and hyperthermophilic bacteria using elastic incoherent 
neutron scattering on INI 3. The optimal growth temp- 
eratures of the different bacteria are 4°C for the 
psychrophile, 37°C for the mesophiles, 75 °C for the 
thermophile and 85°C for the hyperthermophile [21]. 
Interestingly, the average macromolecular flexibility 
was found to be very similar at the opti mal growth 
temperatures of the bacteria with a \/(u 2 ) value close 
to 1.2 A in support of an equivalent state hypothesis: 
i.e. a selection by evolution of similar average dynamics 
for the protein at physiological condition [21]. The resi- 
lience or 'stiffness' of the macro molecules increased from 
the psychrophilic to the hyperthermophilic macromol- 
ecules. The obtained effective force constants ranged 
from (k') = 0.22 + 0.03 N m -1 for the psychrophilic 
Aquaspirillum arcticum, to 0.43 + 0.01 N m -1 and 
0.40 + 0.01 Nm -1 for the mesophilic Escherichia coli 
and Proteus mirabilis and to 0.68 + 0.11 N m -1 and 
0.61 + 0.01 N m -1 for the thermophilic Thermus ther- 
mophilus and hyperthermophilic Aquifex pyrofilus cells 
[21]. In a following study by the same authors on 
rabbit lactate dehydrogenase ( LDH ) and a hyperthermo- 
philic archeal homologue, a \J (u 2 ) value of around 1.5 A 
was reported at the corresponding temperatures of 
optimal catalytic activity. The rabbit LDH and the 



hyperthermophilic archeal homologue were characterized 
by a soft force constant of (kf) = 0.15 N m -1 and a very 
stiff force constant of 1.5 N m -1 , respectively [22]. 

We observed similar equivalent state behaviour for 
vertebrate Hb in our study with an adaption to physio- 
logical temperature. Platypus Hb has the softest 
structure, human Hb is stiffer and chicken Hb was 
found to be the stiffest protein. Crocodile Hb was 
found to be comparatively stiff, with a force constant 
between human and chicken Hb (compare table 2). 
The resilience of the Hb proteins is closer to that of 
rabbit LDH than to the resilience of the archeal or 
the bacterial proteins [21,22]. Stronger internal forces 
are indeed correlated with higher melting temperatures 
of the Hb of all species, and with the average body 
temperatures of Hb of the endotherms, shown in 
figure 4. Platypus Hb, which has the smallest effective 
force constant, shows the larges MSD of all measured 
proteins. However, both crocodile and chicken Hb, 
which are the stiffest proteins in this study, have (u 2 ) 
values that are above those of the softer human Hb. 
The observations confirm that resilience and MSD are 
independent parameters related to the particular 
shape of the force field and that a general assumption, 
for example, that a stiffer protein is always less flexible 
than a softer protein is not valid. 

We now compare the measured RMSD of the Hb 
molecules at the reported body temperatures T h of the 
animals (table 2). The RMSD of human and chicken 
Hb at the respective body temperatures are very close 
to 1.2 A. The RMSD of platypus Hb at T h is slightly 
larger, 1.4 A with a large error. These results suggest 
that a RMSD value close to 1.2 A at T h seems to be 
the optimal flexibility of the proteins at that tempera- 
ture. RMSD of crocodile Hb varies between 1.2 A at 
25°C and 1.3 A at 34°C, which is remarkably close to 
the average flexibility of 1.2 A at the body temperatures 
of the other animals. Strong forces within crocodile Hb 
might keep the amplitudes of motion close to the opti- 
mal flexibility of the protein over the whole body 
temperature range of the animal. The obtained RMSD 
of Hb at body temperature is nearly identical to the 
RMSD of the macromolecules in the bacteria [21], but 
slightly smaller than the RMSD of the rabbit LDH 
and the hyperthermophilic archeal homologue [22]. 

To investigate the role of hydration, we studied protein 
fluctuations of human and chicken Hb at different 
hydration levels. Hydrated protein powders at a hydration 
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Figure 5. Mean square displacements (u 2 ) of the hydrated Hb 
powder samples as a function of temperature. The dotted, 
dashed and solid lines are linear fits to the data in the temp- 
erature range from 10 to 100 K, from 100 to 200 K and from 
252 to 292 K, respectively. The straight lines serve to guide 
the eye. Triangles, human Hb; diamonds, chicken Hb. 



level of 0.4 g D 2 0/1 g dry protein corresponding to 
around one hydration layer are traditionally used to 
study protein dynamics using incoherent neutron scatter- 
ing [41,42] . One hydration layer is assumed to be sufficient 
to activate anharmonic motions above around 180- 
240 K, while global protein diffusion is absent, which 
facilitates the interpretation of the measured MSD. The 
powder samples were measured using the cold neutron 
spectrometer SPHERES. The instrument has a high- 
energy resolution corresponding to observable motions 
in a space- time window of a few angstroms in a few nano- 
seconds. Representative measured data of both samples 
are given in electronic supplementary material, figure 
S2. The determined MSD of human and chicken Hb are 
shown in figure 5. Three linear regions are visible in the 
(u 2 ) as a function of temperature. The first inflection at 
around 100 K is caused by methyl group rotations 
which enter the resolution of the spectrometer [43,44]. 
The second inflection starting at around 200 K is related 
to the so-called dynamical transition. The dynamical 
transition of proteins was attributed to the onset of anhar- 
monicity at 200-240 K, which is related to the melting of 
glassy water molecules in the hydration shell [41,45]. It 
was pointed out that the dynamical transition can also 
be explained by a finite-resolution effect of the neutron 
spectrometer, when the observable protein motions 
enter the time window of the instrument [46-48]. In our 
work, we calculated effective force constants (k') from 
the linear slopes of the (u 2 ) versus temperature in the 
three linear temperature ranges (table 2). Within the 
error bars, the (u 2 ) at all temperatures and the force con- 
stants of human and chicken Hb as hydrated powder 
samples are identical. The conclusion can only be that 
one hydration layer is not enough to activate all motions 
of Hb in the nanosecond time scale which might be rel- 
evant for biological function. The observed difference in 
the dynamics between human and chicken Hb above 
270 K is only activated at a hydration level larger than 
one hydration layer. 
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Figure 6. (a) Force-constant profile for the tetrameric biologi- 
cal unit of chicken Hb (black line), the black squares denote 
the main rigidity peaks in the first subunits and correspond 
to the following residues: Met32-B13, Tyr36-B17, Val96-G3, 
Lys99-G6 and Aspl26-H9 for al; Tyr35-B17, Trp37-C2, 
Glul01-G3, Argl04-G6 and Glnl31-H9 for pi. Force-constant 
variation Ak= k hnuian H b &chicken o Hb of the biological unit 
(red line). Note that 1 kcalmol" 1 A" 2 = 0.7 N m" 1 . (b) Car- 
toon representation of the biological unit of chicken Hb, 
with the al, pi, a2 and P2 units plotted in red, blue, grey 
and orange, respectively. Green van der Waals envelopes cor- 
respond to the central rigid core (formed by G3 and G6 
residues from all four units) and the remaining envelopes 
denote the position of rigid residues at the al/pl and a2/ 
P2 interfaces. 



3.4- Coarse-grained Brownian 
dynamics simulations 

To obtain a detailed understanding of the experimental 
observations on a residue level, supporting coarse- 
grained Brownian dynamics simulations of chicken 
and human Hb were performed. The structures of platy- 
pus and crocodile Hb are not yet available. 

The rigidity profile for the complete tetrameric unit 
of chicken deoxy Hb is represented in figure 6 a (black 
line). As could be expected for this multi-domain 
system, the periodicity of the profile reflects the fourfold 
symmetry in the protein's structure and the profile is 
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Figure 7. (a-d) Force-constant profiles for the four monomeric units of chicken Hb (black lines) and force-constant variations 
A&= Airman Hb — thicken Hb fc> r the monomers (red lines). The red secondary structures plot at the top of (a, b) indicate the 
localization of a-helices along the sequences of the a- and p-units (1 kcal mol -1 A -2 = 0.7 N m -1 ). 



dominated by rigid residues located at the interdomain 
regions [13,49]. The main rigidity peaks correspond to 
residues Val96-G3, Lys99-G6 from the al/2 units and 
residues Glul01-G3 and Argl04-G6 from the pl/2 
units, which form a tight interacting network in the 
core of the biological unit (figure 66, residues in 
green). Secondary rigid areas correspond to residues 
located around Met32-B13, Tyr36-B17 and Aspl26- 
H9 in the al/2 units and Tyr35-B17, Trp37-C2 and 
Glnl31-H9 from the (31/2 units, which interact at the 
al/(31 and a2/(32 interfaces (figure 66, residues in 
red, blue, grey and orange). The variation of the 
force-constant profile when moving from chicken to 
human deoxy Hb (red line in figure 6 a) is mainly due 
to a strong decrease (over —50 kcal mol 1 A 2 = 
35 Nm -1 ) in the rigidity of residues Glul01-G3 and 
Argl04-G6 from the (31/2 units, thus suggesting that 
human Hb is more flexible than its chicken counterpart, 
notwithstanding the fact that human Hb presents a 
slightly denser structure than chicken Hb, reflected in 
the reduced models constructed for both systems, 
which are made of 19 930 and 19 402 springs for 
human and chicken Hb, respectively. We can also 
note that all these residues are conserved in the chicken 
and human Hb sequences, with the exception of Tyr36- 
B17 from the a-units of chicken Hb, which is replaced 
by a Phe residue in the human sequence. Since Tyr 
and Phe sidechains have similar sizes, this mutation 
cannot, however, be detected in our coarse- grain 
model for proteins. 

Using the domain separation procedure detailed in 
Sacquin-Mora & Lavery [13] and Sacquin-Mora et al. 
[27], it is also possible to obtain a rigidity profile for 



each of the four monomers composing the tetrameric 
biological unit. The resulting profiles are shown in 
figure 7 and are qualitatively similar to those obtained 
in the generic study on globins performed by Bocahut 
et al. [12]. The analogous aspect of the profiles for the 
a- and (3-subunits reflects the typical a-helical globin 
fold, with a-helices appearing as grouped rigidity 
peaks along the protein sequences and flexible regions 
in-between denoting, in particular, the CD and EF 
loops. With the exception of Leu66-E15 in the a2 
subunit, all the noticeable mechanical variations 
that can be observed when moving from chicken to 
human Hb correspond to residues from the p-subunits 
becoming more flexible, such as Gly24-B6, Ala27- 
B9, Gly64-E8, Val67-Ell, Phe71-E15, Leul06-G8, 
Glyl07-G9 and Leull0-G12. Interestingly, all the 
aforementioned residues (which are conserved in 
the chicken and human Hb sequences) were previously 
identified as 'Mechanically Sensitive' residues that 
are located at the frontiers between internal cavities 
of the globin fold [12]. In particular, all four members of 
the globin mechanical nucleus (corresponding to pos- 
ition Ell, E15, G8 and G12) present an important 
decrease of their rigidity in the (31 unit of human Hb 
compared with chicken Hb. 

Altogether, the mechanical variations when moving 
from chicken to human Hb are mainly located at the 
central interface of the biological unit (figure 66), on 
mechanically sensitive residues from the monomers, 
and only concern residues that are conserved in both 
sequences. Therefore, the thermal adaption of Hb 
internal dynamics might result from a long-range muta- 
tional effect, similar to what was previously observed in 
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the bacterial photosynthetic reaction centre [14], where 
point mutations would affect the mechanical properties 
of residues located more than 15 A away. 
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